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U . ABSTRACT 

^, ■ We explore the design of a space mission called Project Lyman that has the goal of quantifying the ionization history of the 

universe from the present epoch to a redshift of z ~ 3. Observations from WMAP and SDSS show that before a redshift of 
z > 6 the first collapsed objects, possibly dwarf galaxies, emitted Lyman continuum (LyC) radiation shortward of 912 A 
K^ I that reionized most of the universe. Theoretical estimates of the LyC escape fraction (fesc) required from these objects to 

complete reionization is fesc ^ 10%. How LyC escapes from galactic environments, whether it induces positive or negative 
feedback on the local and global collapse of structures, and the role played by clumping, molecules, metallicity and dust are 
p^ ■ major unanswered theoretical questions, requiring observational constraint. Numerous intervening Lyman limit systems 

I ^ I frustrate the detection of LyC from high z objects. They thin below z '--^ 3 where there are reportedly a few cases of 

("^ ■ apparently very high fesc- At low z there are only controversial detections and a handful of upper limits. A wide-field 

OO ! multi-object spectroscopic survey with moderate spectral and spatial resolution can quantify fesc within diverse spatially 

^D ■ resolved galactic environments over redshifts with significant evolution in galaxy assemblage and quasar activity. It can also 

calibrate LyC escape against Lya escape, providing an essential tool to JWST for probing the beginnings of reionization. 
. , ■ We present calculations showing the evolution of the characteristic apparent magnitude of star-forming galaxy luminosity 

rS I functions at 900 A, as a function of redshift and assumed escape fraction. These calculations allow us to determine the 

jrt ■ required aperture for detecting LyC and conduct trade studies to guide technology choices and balance science return 

against mission cost. Finally we review our efforts to build a pathfinding dual order multi-object spectro/telescope with a 
(0.5°)^ field-of-view, using a GSFC microshutter array, and crossed delay-line micro-channel plate detector. 
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1. INTRODUCTION 

Observations show that by a redshift of z « 6 most of the universe passed through an era of reionization precipitated by an 
increase in the metagalactic ionizing background (MIB).^ Star-forming galaxies and quasars are generally acknowledged 
as the most likely sources responsible for reionization, but how their relative contributions change over cosmic time is the 
subject of considerable debate.^ 

Number counts of galaxies and quasars at z '^ 6 show that low mass star-forming galaxies vastly outnumber the quasar 
population, ''^'^ strongly suggesting but not proving, that star- forming galaxies dominate reionization.^ The faint end 
slope of the quasar luminosity function near z ~ 6 is poorly constrained, so it is not clear if quasars can produce ionizing 
radiation on par with that required from galaxies to sustain reionization epoch. ^ There are also indications that the initiation 
of reionization may require a "hard" spectral energy distribution (SED) or an fe > 0.2 that increases with redshift.^' ^ 
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The relative contributions to the MIB for each population depends upon the faint end slope of the luminosity function, 
the escape efficiency fesc of ionizing radiation emitted shortward of the Lyman edge at 91 1.7 A, the H II clumping factor 
C =< p^ > / < p >^ and the spectral hardness. ^'^ The evolution of these quantities has a fundamental influence on the 
state of the MIB, the intergalactic medium (IGM) and the formation of structure throughout cosmic time.'" 

Ionizing radiation can produce either positive or negative feedback on the star formation process on global and local 
scales. The nature of the feedback depends in part on the clumping environment. '^'^^ Positive feedback occurs when 
photo-electrons catalyze the formation of H2, which allows the gas to efficiently cool and collapse. Negative feedback 
occurs when photoionization heating temporarily halts gas collapse by increasing the Jeans mass. How feedback relates 
to the conditions that favor the escape of Lyman continuum (LyC) photons from star-forming environments and the role 
played by clumping is not well understood. We know escape of LyC photons from wavelengths shorter than the H I 
ionization edge at 91 1.7 A is very difficult,'"' '^ some have even said impossible.'"* Mean H I column densities for normal 
galaxies are Nhi = 10^' cm^^, yet it takes only a column of 1.6 x lO'^ cm^^ to produce an optical depth of unity at the 
Lyman edge. 

For quasars this does not seem to be much of a problem. Below the Lyman edge quasars emit a "hard" SED, which 
can ionize both H I and He II, while star-forming galaxies, which emit a "soft" SED, are capable of only ionizing H I. The 
characteristic ultraviolet luminosity of quasars is 10 times that of the galaxies, so the IGM surrounding quasars tends to be 
devoid of neutral material. There are numerous examples of quasar SEDs that show no sign of a break in the continuum 
emission through the Lyman edge.'^''^ Similar observations for galaxies are rare'^"'^ and often controversial. ^""^^ It is 
fair to say that there is considerable uncertainty regarding the physical conditions that favor LyC from star-forming regions 
at any epoch. 

It is somewhat simplistic to characterize reionization as a process caused by quasars with /^ = 1 or star-forming galaxies 
with fe < i- The central engines of active galactic nuclei (AGN) come in a variety of masses, have intermittent duty cycles 
and can be obscured by host galaxy environments, resulting in a softer SED.^^ Some fraction of quasars do exhibit a break 
at the Lyman edge,^^ contrary to the conventional assumption of unity escape fraction.^ Furthermore, it is now accepted 
that black holes reside in the nuclei of most if not all quiescent galaxies, ^^'^^ so the effects of previous AGN activity within 
an apparently dormant environment must also be considered. 

Observations showing the transition of black H I absorption troughs to a forest of narrow absorption features at redshift 
z <; 6,' compel us to infer that LyC photons somehow escape star- forming environments and hence must be detectable 
at some level. A failure to observationally confirm LyC escape at any level would necessitate a serious reevaluation of 
conventional theories for the reionization era. Unfortunately, direct measurements of the LyC becomes difficult for z J> 3 
due to the increasing optical depth of the IGM.'^ The collection of Lya forest systems (N(H I) <; lO'^ cm^^), Lyman limit 
systems (lO'^ <; N(H I) <^ lO'^ cm^^) and damped Lya systems (lO'^ cm^^ <^ N(H I)) all increase toward high z- The 
probably that a galaxy will be obscured by an intervening system as a function of redshift has recently been quantified by 
Inoue & Iwata.^** They characterize the detection probability in terms of the magnitude difference between the continuum 
redward of the Lyman edge and the detection limit blueward of the edge. For a 2.5 magnitude difference (a factor of 10 in 
flux) they find detection probabilities of 98, 92, 73, 19, 0% at z = 1, 2, 3, 4, 5. From this we conclude that stastical studies 
to directly examine the physics of LyC leak are best conducted at z ^ 3. Detectable LyC leakers beyond z ^ 4 will be rare. 

Here we discuss the requirements for detecting LyC leakage and explore the parameter space for a spacebased survey 
mission called Project Lyman, utilizing a wide-field spectral/imaging technique. The goal of this project is to quantify 
the ionization history of the universe at redshifts z <^ 3 by determining the fg^c as a function of galactic morphology and 
clustering along with the local environmental factors of clumping, metallicity and gas-to-dust ratio. We will search for low 
redshift analogs to the low mass star-forming galaxies thought to be responsible for initiating reionization, and determine 
whether Lya emission can serve as a proxy for the LyC escape from high redshift galaxies. 

We first briefly review the star formation history of the universe, discuss the expectations for how the escape fraction 
evolves, and give a plausibility argument for the existence of low redshift analogs to high redshift reionization sources. A 
review of the fesc detections and lower limits will be presented. We also discuss how James Webb Space Telescope (JWST) 
reionization studies could benefit from the calibration of Lya escape with LyC escape that this project proposes to make. 
We then use the redshift luminosity functions presented by Arnouts^^ to estimate the areal densities of objects and the 
required detection sensitivities as a function of redshift. The aperture requirement for a long duration spacebased mission 
is determined and compared to the capabilities of a prototype spectro/telescope we are fabricating for a series of sounding 
rocket missions to provide proof-of-concept for the optical design. 



2. STAR FORMATION HISTORY OF THE UNIVERSE 

As the universe evolves the dark matter halos, which make up the bulk of the matter in the universe, start to coalesce and 
attract baryons that turn into star filled luminous galaxies. Early on, the mass of these halos and the baryon number within 
them is small and the star production is low. As the merger process accelerates, the characteristic absolute luminosity in 
the ultraviolet rest frame, a proxy for the star formation rate, rises from M*^, « -20 at z « 6-7 and reaches a peak of M*^ « 
-21 after about 2 Gyrs near a redshift of w 3.5. In the ensuing 11 .7 Gyrs M*j, has fallen to « -18.^ The intensity of the MIB 
closely follows the stellar luminosity density. 

Merger studies suggest that the peak luminosity near z ^ 3 is associated with the assemblage of the most massive 
galaxies, undergoing an enormous burst of star formation. This peak also roughly coincides with the peak in the quasar 
number density, and maybe associated with the He II reionization era near z ^ 3.^'^ These most massive galaxies quickly 
form stars and perhaps exhaust their fuel reservoirs in the process. Today they are dominated by an older low mass 
population of stars and are members of the "red sequence" in the color bimodality distribution of galaxies. ^"^ The rapid 
star formation period of these most massive galaxies is somewhat analogous to the short nuclear burning timescales of the 
most massive stars. 

It has been observed that in the low redshift universe following the star formation peak at z ^ 3, activity has since 
shifted from high mass systems to low mass systems in a process known as "downsizing."^''^^ In other words most of 
the present day star formation is occurring in the lower mass "blue sequence" galaxies. It is possible that we are now in 
an epoch where the lowest mass halos, which at a z '^ 7 were not clumpy enough to support star formation, have gained 
enough mass to initiate stellar collapse. The chemical enrichment in the low mass halos at low z is an open question. It 
depends on the efficiency of chemical feedback from galaxies into the the IGM. It maybe that the low metallicity, low 
mass star burst systems of today are analogs to the numerous low mass, low metallicity dwarf galaxies thought to have 
been responsible for sustaining and possibly initiating the reionization epoch. Detecting LyC escaping from the most metal 
poor, low luminosity systems in the low redshift universe will tell us much about escape process in the early universe. 
These objects are expected to be rare, on the faint end slope of the luminosity function, and will require wide field survey 
techniques to detect. On the other hand, they are expected to have higher escape fractions on average, which will aid in 
their detection. 

3. ESCAPE FRACTION EVOLUTION - EXPECTATIONS AND OBSERVATIONS 

In the high redshift universe (z ~ 6) low mass star-forming galaxies are expected to vastly outnumber the quasar popula- 
tion,'"* yielding a soft MIB and an early H I reionization era, as long as the average the LyC escape fraction is fesc ^ 10 - 
20%.'*'^'^^ There are no direct observations of LyC escape in this epoch. 

Towards intermediate redshifts (z ^ 3) an increasing quasar population is expected to harden the MIB and usher in an 
era of He II reionization. However, comparative observations of the H I and He II Lya forests over the redshift range 2.3 

< z < 2.9 show that in this era star-forming galaxies still supply a significant fraction of the MIB photons.^'*' ^^'^^ This 
conclusion is supported by the Steidel et al. '^ detection of LyC escape in a composite spectrum of 29 Lyman break galaxies 
with a mean redshift of z = 3.4, suggesting an escape fraction significantly higher than 5%. More significantly, Shapley et 
al.'^ have found 2 Lyman break galaxies (LBG) at intermediate redshift (z ^ 3) with remarkably high escape fractions 40 

< fesc < 100%. There is also a new report of LyC detected from 10 Lya emitters (LAEs) and 6 LBGs near z = 3. 1, using a 
narrow band filter technique.'^ The average fe > 15% estimated for this sample is somewhat tentative, pending follow-up 
spectroscopy on a few of the LAEs. 

In the low redshift era, calculations indicate that the contribution of star-forming galaxies to the MIB will continue to be 
competitive with that from quasars if on average fesc ^5%.^ Thus if star- forming galaxies are indeed major contributors to 
the MIB, then a non-trivial fraction of LyC photons should be escaping and directly detectable with spectroscopy, provided 
the source is at a redshift high enough to escape the H I absorption "shadow" of the Milky Way; z ^ 0.02. The LyC 
absorption shape should have an envelope <=<: X^^ shortward of the Lyman edge.^^ Low redshift observations have been 
less successful than those at z ^ 3 in detecting LyC. A recent low redshift search near z ^ 1.3 using far-UV images of the 
Hubble Deep Field-North and the Hubble Ultra Deep Field reports fe < 0.08, using a stack of 21 candidates.'^ 

Leithereret al.^" found upper limits on fesc < 1-15% from Hopkins Ultraviolet Telescope measurements of 4 starburst 
galaxies with redshifts 0.018 < z < 0.03. Hurwitz et al.^' reexamined the same data and found upper limits fesc < 3 - 57%. 



Deharveng et al.^^ used the Far Ultraviolet Spectroscopic Explorer (FUSE) to observe Mrk 54 at z = 0.0448 and found 
fesc < 6%. Bergvall et al.^^ also used FUSE to observe the blue compact galaxy (BCG) Haro-11 at z = 0.021. Their 
detection of 4 < fesc < 10% was questioned by Grimes et al.^^ who found only an upper limit fesc ^ 2% in a reanalysis of 
the data. 

We conclude that a meaningful low z program needs to explore the regime fesc < 5%, with a sensitivity well below the 
ft/5£ background limit ~ 10^'^ergscm^^s^' A^' = \FEFU (femto-erg flux unit). 

4. ENVIRONMENTAL INFLUENCES ON LYC AND Lya ESCAPE RELATIONSHIPS 

The escape fraction at z '-^ 3 appears to have been higher than at low z. Perhaps this is due to a paucity of well formed 
neutral disks around galaxies in the earlier epochs. We know the Hubble sequence began to form at z ~ 1 - 1.5^^ when the 
expansion of the universe caused the rate of major mergers to slow, allowing for the first time massive H I disks to settle and 
remain stable on global scales. Galaxies at z — 2-3 are likely to have high escape fractions, because they have not had time 
to form stable H I disks, due to the combined effects of a much higher major merger rate^^ and starburst- or AGN-driven 
outflows.'^ The combination of tidal disruption and ionization of the residual gas in young star clusters, perhaps aided 
by a strong MIB, makes it easier for the LyC emission from star formation to exceed the recombination rate^** and create 
"density bounded" H II regions with high escape fractions. The low densities within the disrupted environments also aid 
large scale stellar and supernovae generated outflows to create superbubbles and chimneys through which LyC (and Lya) 
can escape."*' 

This will be much less the case at lower redshifts (0 < z < 1 .5), where giant H I disks have begun to settle and outflows 
are significantly smaller Lower fesc may result, due to the relatively high density of H I surrounding "radiation bounded" 
H II regions. Nevertheless, we expect a trend where metal poor dwarfs and irregulars have higher fesc than bulge and disk 
type galaxies because they tend to reside on the outskirts of galactic clusters and are surrounded by a more tenuous IGM 
with a higher ionization fraction on average. These objects may be the low z analogs to the low mass LyC emitting objects 
at high z. Furthermore, examples abound at low z of mergers and objects with unstable disks, large scale outflows, and 
young star clusters embedded within clumpy multi-phased media. In short, all the same physical conditions conducive 
to LyC escape at high z can be studied directly in low z galaxies, up-close and in-depth, with perhaps a somewhat higher 
metallicity. Nevertheless, by characterizing fesc at low redshift in a diverse sample of galaxies we gain insight to the physics 
of LyC escape, with which we can constrain the contribution of high z galaxies to the ionization history of the universe. 

4.1 Is Lya Escape a Proxy for LyC Escape? 

Spectroscopy provides the most direct means to measure LyC and Lya over a continuous range of redshifts and presents a 
unique opportunity to determine whether a proxy relationship exists between Lya and LyC escape. This effort is extremely 
important to the James Webb Space Telescope {JWST) key goal of identifying the sources responsible for initiating and 
completing the epoch of reionization. Moreover, Lya emission is thought to be a beacon for the formation of structure in 
all epochs.^^ It is the primary ionization diagnostic available to JWST and may provide a detailed view of the beginning of 
the reionization era at redshifts z > 6.^^ 

At first glance, it appears that the probability for Lya photons to escape is even lower than for LyC photons. The line 
core optical depth reaches unity for a column of N{HI) = 1.4 x lO'^ cm^^ with a doppler velocity of 10 km s^'. Case B 
recombination theory posits that Lya photons created in a recombination event are reabsorbed "on-the-spot," and become 
trapped in "radiation bounded" H II regions. In the absence of dust they scatter endlessly in successive absorptions and 
re-emissions until they random walk into the line wings and escape. In the process a characteristic double-peaked profile 
is created with a self-absorbed line core.'*'* Adding dust further blackens the core and reduces the peaks. In contrast, a 
single absorption of a LyC photon results in its destruction, and the creation of a Lya photon. Adding velocity fields 
further aids Lya escape by doppler shifting the photons out of the optically thick line core and into the thin line wings, 
a process which is unavailable to LyC photons. Outflows skew double peaked profiles towards the red, inflows towards 
the blue.^^'"*^ However, velocity gradients can aid LyC escape by creating low density regions with lower probability for 
continuum absorption. 

Scattering of Lya through a multi-phase (clumpy) medium has been found to enhance Lya equivalent widths (EW)."*^ 
Examples of Lya EW boosting may have been observed in high z star-forming systems, which often exhibit impressively 



high EW > 150 A.^*' A more typical value is EW « 20 A (rest frame) as found in the composite LBG spectrum at z « 3.4,'^ 
which is similar to that found at low z in star-forming galaxies. ^^ 

It is not clear whether these processes will lead to a correlation or an anti-correlation of LyC and Lya escape. We note 
the intensity of Lya is proportional to the number of LyC photons that do NOT escape. A simple formula, following Mao 
et al.,^^'^*' for calculating the luminosity of escaping Lya photons is, 

LLy, = \Q{^-fesc)fLy, (D 

where fiy^ is the escape fraction of Lya photons and Q is the total LyC emission rate. Ionizing radiation (Q) has to be 
present to create Lya in the first place, but if //,,„ is roughly constant and is mostly independent of the physical processes 
that cause LyC leakage, then as fesc increases L^y^ decreases. On the other hand, if the Lya and LyC escape processes are 
co-dependent (i.e. fesc ~ fhya = /) then at / w 0.5 we find a peak luminosity for L^ya = \Q- 

Commensurate observations of Lya and LyC to date are suggestive but not clearcut. In the Shapley et al.'^ data on 
14 star-forming galaxies at z ^ 3 there are 5 objects without LyC and Lya, 7 objects with Lya and no LyC, 1 object with 
LyC {fesc ^ 50%) and strong Lya, and 1 object with LyC {fesc > 65%) and marginal Lya. Lya and LyC emission do 
appear together, but there is some variation. Further efforts to characterize this relationship will require a large statistical 
sample. A quantitative assessment of the contribution of galaxies to the ionizing background will require a large number 
of observations before LyC and Lya luminosity functions can be established.^' 

4.2 Project Balmer 

In support of this effort. Project Balmer, a groundbased observing program to acquire spectra of Ha and HP, will provide 
a measure of the dust attenuation at visual wavelengths. It will also allow the determination of the absolute luminosity 
of the intrinsic LyC produced within, but not necessarily escaping the star-forming galaxies. After accounting for the 
effects of dust attenuation and temperature, using standard emission line diagnositics,^^ the Lya escape fraction can be 
found by comparing the observed Lya/Ha flux ratio to the intrinsic recombination flux ratio. Together Project Lyman and 
Balmer provide a complete accounting of the LyC and Lya escape budgets for addressing the LyC escape and Lya escape 
correlation. If such a correlation can be established with low redshift observations, then the relatively easy to observe Lya 
could be used to quantify the contribution of star-forming galaxies to the MIB at all redshifts. 

5. LYC DETECTION REQUIREMENTS ^ Z ^ 3 

GAZ^Xhas shown that there are hundreds of far-UV emitting galaxies in each square degree of the sky. We have suggested 
that a wide field spectroscopy survey can efficiently search for LyC leak.^^ To plan a search we need to determine the 
required instrumental sensitivity to detect Lyman continuum escaping star-forming galaxies representing a diversity of 
morphological types over the redshift interval < z <; 3. The top level instrument design will flow down from the required 
detection sensitivity. Since these observations must be made in space it is important to have a realistic estimate of the 
scientific potential as a function of aperture to control mission cost. To carry out this assessment we need to know the areal 
density of Lyman continuum candidates on the sky, as a function of apparent magnitude (i.e. observer frame flux) and 
redshift. 

In Figure[T]on the left we plot an estimate based on luminosity functions provided by Arnouts et al. using the Schechter 
parameterization,^^ for the redshift intervals indicated in the caption. The figure shows the logarithm of the number of 
galaxies per unit magnitude per square degree as a function of the restframe ultraviolet {^ 1500 A) AB magnitude as 
measured in the observer's frame. We converted from the number per comoving volume to the number per square degree 
by multiplying by the comoving volume per solid angle, calculated for the mid-point in each redshift interval. Each curve 
has an * to indicate the location of the characteristic magnitude 'Wisogn , ,) where the Schechter function makes the transition 
from the exponential cutoff in galaxy counts at the bright end of the luminosity function to the power law extension of the 
faint end. The interval for the ordinate of each curve spans 5 magnitudes (a factor of 100 in flux) centered on 'Mi5oo(i+,)- 

We used the formula^^ 
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Figure 1. LEFT PANEL - Surface densities as a function of observer's frame apparent magnitude for galaxy populations with redshifts 
between - 0.2, 0.2 - 0.4, 0.4 - 0.6, 0.6 - 0.8, 0.8 - 1.2, 1.8 - 2.3, 2.4 - 3.4, estimated following Arnouts.^^ RIGHT PANEL - 
Apparent magnitude of the LyC emitted by ra*9oo( l+z) galaxies as a function of redshift and assumed escape fractions of /f,sc=0.02, 0.04, 
0.08, 0.16. Strawman far- and near-UV background limits as a function of wavelength for spectro/telescopes with 0.5 m, 1 m, and 2 m 
apertures are overplotted. For reference flux densities (/j^) in FEFU = 10^'^ (ergscm^^ s^' A^') are indicated with long-dashed lines. 
The redshifts are marked at the top of the figure. The far-UV assumes SiC mirrors and a Csl photocathode. The near-UV assumes MgF2 
over Al mirrors and GaN photocathode. 



to transform from absolute magnitude in the 1500 A restframe to the apparent magnitude in the observer's frame, where 
dl{z) is the luminosity distance in Mpc, z is the redshift, and M^y is the absolute restframe UV magnitude. We have 
ignored the K-correction, because it is small in comparison to the other terms in Eq. |2] 

We convert the surface density estimates for the 1500 A restframe magnitude to Lyman continuum magnitudes, using 
the scale factors given by Leitherer^^ for starbursts. For a continuous starburst a typical ratio is /1500//900 ~ 2, assuming a 
solar metallicity and a Salpeter initial mass function with an upper mass cutoff of 100 Mg. This ratio is relatively insensitive 
to age. The range is 1 .5 ^ /1500//900 <; 3 for ages 10 - 900 Gyr. A factor of 2 amounts to Sm^Qg" = 2.5 log/1500 //900 ~ 0.75. 

In Table[T|we list the apparent magnitude of the LyC as a function z and f^sc. 
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where hrtiesc = 2.51og/est-. We list the redshifted wavelength of 912 A at the bottom of the table. Table [T]and Figure[T] 
show that there are hundreds to thousands of galaxies at the 'WgQon 1 \ per magnitude per square degree at all redshift ranges. 
There will be no lack of LyC candidates provided we have enough sensitivity to reach the wLq. j , ■. galaxies. 
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6. INSTRUMENT BACKGROUND LIMITS 

We now turn our attention to determining the background limits imposed by zodiacal light, airglow and detector dark 
counts. We will assume current state-of-the-art efficiencies for the mirrors, gratings and detectors. We will also assume an 
optical design with the absolute minimum number of optical components - one primary mirror, one secondary reflection 
grating, one detector. We will examine two "strawman" telescope configurations to cover the widest available redshift 
range. One is for the far-UV with a bandpass of ^ 900 - 1800 A, incorporating SiC coated optics and a semi-transparent 
Csl photocathode for work at the lowest redshifts (z <; 0.5 for both Lya and LyC, z ^ 1 for LyC). Another is for the near- 
UV with a bandpass of ^ 1800 - 3600 A, using MgF2 over Al optics and semi-transparent GaN photocathode,^^ for work 
at the highest redshifts (z > 0.5). The method for performing the calculations has been given previously.^^ As a starting 
point we will assume the spectral and spatial resolutions with the parameters listed in Table |2] taken from the FORTIS 
instrument, and simply scale the results with the telescope clear area. As the near-UV channel is an "octave" lower than 
the far-UV channel this amounts to using a grating with half the ruling density. 

In the far-UV the detector dark counts and the geo-coronal Lya line are the main contributors to the background. Other 
Lyman series, O I, N I and N? also contribute, but at lower levels. In the near-UV longward of 2000 A the zodiacal light 
becomes progressively more intense. In the Right Panel of Figure [T]we plot 'w^ood +-)(/"') selected from from Table [T] 
The background limits for the far- and near-UV strawmen are overplotted for 0.5 m, 1 m, and 2 m apertures. 

We see that in the near-UV the zodiacal light ceases to be a problem shortward of 2500A. We conclude that aim 
far-UV telescope would be most useful in the 0< z < 0.2 region. It could reach ?m|qqq galaxies with fesc > 0.025. A 2 m 
telescope in the near-UV with a GaN photocathode could reach OTjgQ,j (z « 1) with f^sc ^ 0.03. If the faintest blue galaxies 
have fesc ^ 0.08, then even a 0.5 m telescope could detect galaxies on the faint end of the luminosity function above m*^^^^. 
Detection of LyC leakage for z > 2 will require a telescope in excess of 4 m diameters, to beat down the "rising sun." 
Project Lyman for z > 2 might best be performed from the ground, but additional penalties of atmospheric background and 
attenuation will require even larger apertures. A balloon flying at 40 km might avoid such penalties. 

7. FORTIS PATHFINDER INSTRUMENT - PROGRESS AND PERFORMANCE EXPECTATIONS 

The compelling challenge of detecting LyC leakage in fields of star-forming galaxies provides the impetus to develop a 
highly efficient, multi-object spectrograph for surveying extended regions in the far- and near-UV. Here we briefly describe 
some design aspects, operational considerations and progress towards fabricating a pathfinding instrument that will be 
launched on a sounding rocket to provide scientific and technical proof of the FORTIS (Far-ultraviolet Off Rowland- circle 
Telescope for Imaging and Spectroscopy) instrument concept. A schematic is shown in the Left Panel of Figure|2l 

7.1 Optical Layout 

FORTIS is a Gregorian style spectro/telescope with a ruled secondary mirror It satisfies the requirements of high efficiency 
in the far-UV by having only two reflections and dual order spectral channels. ^^"^^ A zero-order image of the slit plane is 
formed on-axis at the secondary focus and the positive and negative spectral orders are dispersed in a plane to either side. 



Table 2. FORTIS Instrument Summary 



Parameter 


Value 


Parameter 


Value 


Prime Focus Platescale 


206 "3 mm-> 


Inverse Dispersion 


20 A mm-' 


Primary Diameter 


500 mm 


Field of View 


1800" X 1800" 


Primary Radius 


2000 mm 


Microshutters in FOV 


43 X 86 


Central Obscuration 


219 mm 


Single Shutter Area 


15.9" X 36.9" 


Secondary Diameter 


165 mm 


Telescope Clear Area 


1600 cm2 


Secondary Radii xz, yz 


520.83, 520.68 mm 


Spectral Bandpass 


900 -1700 A 


Eccentricity xz, yz 


2/3, 0.66679 


Spectral Resolution 


1500 


Secondarylmage Distance 


1562.50 mm 


Spatial Resolution 


l-4"outto±450" 


Zero-order Platescale 


41.25" mm- 1 


Detector Pixel 


0.012 ^m 


System Focal Ratio 


f/10 


Detector Area 


45 mm x 170mm 






Peak Effective Area 


- 50 cm^ (SiC coating) 
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Figure 2. LEFT PANEL - FORTIS sounding rocket optical layout. RIGHT PANEL - FORTIS raytrace. Point spread functions as a 
function of wavelength and vertical off-axis angle calculated with a geometric raytrace. Inset points are wavelength pairs, separated by 
1 A, expanded by ten and sampled with (0.015 mm)^ bins. Off-axis image rotation will be corrected by distortion mapping. 
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Figure 3. Target simulation for M33. The zero-order image of the target is in the center. The rectangular regions are the slits with the 
brightest integrated signal on the row. Only one slit is allowed per row to avoid spectral confusion. The light from each slit is diffracted 
into the positive and negative orders to either side of the zero-order image. A rich spectral field is evident. 

A combination of holographic ruUng and triaxial figuring of the secondary mirror along with an optimum choice of 
platescale, greatly reduces spectral aberrations and field curvature over a wide field-of-view (« 0?5). The image quality 
in the spectral orders is a vast improvement over the prime-focus Rowland circle mount used by FUSE, which only had 
point source spectroscopic capability. Examples of the point spread functions are shown in Figure |2] Right Panel. The 
Gregorian configuration allows placement of a microshutter array at the prime focus. This creates a powerful mutliobject 
spectrograph with a considerable multiplexing advantage. 

7.2 Spectral Simulation 

We show a spectral simulation of the M33 field in Figure [5] We used the GALEX flux calibrated FUV and NUV images 
of M33 to simulate the emission using the (1''5)^ square pixels provided by the GALEX images. This was done to keep 
the image size of the calculation reasonable as the FORTIS design calls for 0"5 pixels. We assume that each GALEX pixel 
has a fi powerlaw with an slope Po derived from the FUV and NUV images. Given the difference between the observed 
slope and an intrinsic slope (Ap = -2.5 - po) we model the atomic and molecular hydrogen absorption, assuming canonical 
gas-to-dust ratios and molecular fractions. The simulated spectra are produced using 0.73 A bin^^ We used the estimated 
FORTIS effective area curve, ^^ to convert fx to total counts for an integration time of 400 s appropriate to a sounding 
rocket flight. A 3 KRayleigh geocoronal source (9x 10^'^ ergs cm^^ s^' arcsec^^), sets the background limit for each 
slit. The effects of astigmatism have been neglected. The counting statistics have been simulated by adding or subtracting 
the square root of the counts per bin multiplied by a random number drawn from a normalized Gaussian distribution. 
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Figure 4. Schematic cross section of a microsiiutter unit cell showing the key elements. Right - 64 x 128 microshutter array with an 
area of (12.8 mm) is shown at the center, mounted to a Si carrier and bonded to a fanout board in this lab version. The FORTIS fanout 
will minimize the vertical height to ~ 38 mm. 

7.3 Microshutter Assembly 

The microshutter arrays currently being developed at Goddard Space Flight Center (GSFC) for the Near Infrared Spec- 
trograph (NIRSpec) on the James Webb Space Telescope (JWST) consist of a two-dimensional array of closely packed 
clear aperture slits, each with an independently selectable shutter.^^ Clear aperture slits are ideally suited for multi-object 
spectroscopy in the far-UV where traditional techniques relying on optical fibers, lenslet arrays, mirrored slits or image 
slicing cannot be used because of the poor efficiency of all optical materials in this bandpass. The following describes the 
microshutters, their operation, the control electronics and the FORTIS zero-order imager interface. 

7.3.1 Microshutter Description 

Microshutter arrays with 64 x 128 shutters were originally developed as production and lifetime test models for the larger 
175 X 384 cyrogenic devices slated for NIRSpec. The arrays share the same shutter unit cell geometry and electrode 
structures. We will use an uncooled 64 x 128 aiTay on FORTIS. 

An array is composed of a thin shutter membrane and light shield mounted to a 100 fxm thick support grid. The support 
grid is an etched silicon wafer with a 1 00 ^m x 200 ^m pitch of rectangular holes. The matching shutter array is machined 
into a silicon nitride membrane 0.5 /jm thick. An individual shutter blade is suspended from the shutter frame by a torsion 
flexure. Reliability of the shutters has been evaluated and it has been shown that fatigue of the hinges is not a credible 
failure mechanism. Tests demonstrate that the number of failed open shutters equals < 1 % after 10^ cycles. Shutter blades 
are surrounded by an aluminum light shield, which overlaps the edges. This provides high on-to-off contrast ratio ^ 7000: 1 
The unit cell open area ratio is 0.69.^" 

The shutters are coated with a layer of high permeability magnetic material to allow for magnetic actuation of the 
shutters. Electrodes, deposited on one of the interior walls of the shutter support grid and the shutters themselves, allow 
for electrostatic latching and release. Figure H] shows the unit cell with exaggerated dimensions, indicating the different 
structural and electromechanical layers. 

An aperture is opened by rotating the shutter by 90° on the torsion flexure. A magnetic actuation and electrostatic 
latching method is used. Shutters are rotated out of the focal plane by the torque created as an external magnet is scanned 
over the shutter array. They are held open by an electrostatic potential applied between the open shutters and electrodes 
deposited on the interior wall of the support grid. To enable the electrostatic latch an aluminum coating is deposited on 
the microshutters themselves to form address rows. Address columns are formed by depositing backside electrodes on the 
walls of the frame in a direction orthogonal to the rows, using angled deposition beam. During actuation and latching of 
the shutters, H-24 V is applied to the shutters, and 24 V is applied to the wall electrode. Either the row voltage or the column 
voltage is sufficient to hold the shutter open. A shutter is released when both the row voltage and the column voltage are 
set to zero. This allows full random access, using a two dimension crosspoint address scheme. 

7.3.2 Microsiiutter Components 

A set of arrays and control electronics will be fabricated by GSFC. Selected aiTays will be "bump bonded" to a Si aiTay 
carrier. The Si array carrier will be mounted to a custom designed printed circuit (PC) fanout board with a small vertical 
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Figure 5. Left Panel - Far-UV zero-order cylindrical achromat. Middle Left - Zero-order raytrace PSFs at the focal plane extremes 
without correction. Middle Right - Zero-order corrected with the achromat. One pixel = 0.012 /jm = Of 5. Right Panel - A rendering of 
the detector assembly showing the position of the zero-order cylindrical achromat. 

extent, w 38 mm, to minimize obscuration of the beam reflected off the secondary. The fanout board provides the paraUel 
electrical connections to the rows and columns of the individual shutters to allow application of the ± 24 V levels for 
opening and closing the shutters from a high-voltage shift register (HV584) fabricated by SuperTex. The HV584 is a 
128 channel serial to parallel converter with push-pull outputs. The HV584 serial inputs come from level shifted logic, 
controlled by 5 volt CMOS inputs. The appropriate serial waveforms are generated by a microprocessor. Two HV584s are 
required, one for the rows and one for the columns. 

7.3.3 Zero-Order Microshutter Interface - ZOMI 

A zero-order microshutter interface, between the HV584 level shifter and the zero-order imaging channel, will allow 
addressing of the individual microshutters during flight. It consists of a an onboard microprocessor to create a bit map 
array of the zero-order target field, a field programable gate aiTay (FPGA) to pass the bit map to the 584, and a magnet 
sweeper for actuating the shutters. FPGAs will be programmed to take a bitmap aiTay and output a serial waveform to the 
584 interface timed with the magnet sweep actuation. The FPGA will also have various functions preprogrammed into it 
such as, open all shutters or close all shutters. The magnet will be mounted on a linear translation stage thats travel back 
and forth, for opening and closing, by « 1.5 inch. 

7.3.4 Zero-order Image Correction 

Astigmatism can be corrected holographically, but in only one order at at time. Postive and negative order astigmatism 
correction, as required by a dual order design, can only be carried out with a grating figure that has slightly different sagittal 
and tangential radii. ^^ This surface creates a point source image with an astigmatic height of « 280 ^m (11 f'5). 

We have designed an achromat using CaF2 and MgF2 cylindrically shaped lenses to correct for astigmatism in the 
imaging focal plane. The use of CaF2 has the added benefit of eliminating geocoronal Lya in from the imaging region. 
In Figure|5]the lens design is shown in the Left Panel, along with the raytrace results showning the imaging improvement 
provided by the achromat (Middle Panels). For this raytrace input point spread functions (PSFs) were radial Gaussian spots 
with 1" FWHM and assumed a flat spectral response from 1300 - 1900 A. Corrected PSFs are « 2"over most of the 1800" 
X 1800" FOV. 

7.4 Detector 

The micro-channel plate (MCP) detector assembly is in the Right - Panel of Figure |5] The anode is a crossed delay line 
with three separate readouts. A 42 x 42 mm^ readout for zero-order, and two 64 x 42 mm^ readouts for the spectral 
regions. There are 4 mm gaps between the anodes. The MCPs also are separated by 4mm gaps but share a common output 
voltage plane. The output plane to anode voltage is maintained constant with zener diode. The input voltages are separate 
for each MCP to allow the spectral channels to be turned off when acquiring the zero-order image. Otherwise geo-coronal 
Lya would swamp the spectral channels with counts, which, while not fatal, would be too high a count rate to be accurate. 

The detector housing has a door, which will open in flight and during preflight testing, to allow the light to reach 
the detectors. The detector assembly will also include, motors to actuate the door, a pump-out port, ion and getter pump 
assembles (not shown in the drawing). The entire spectro/telescope is housed in vacuum tight skins and evacuated prior to 
flight. 



7.5 Target Acquisition 

The (1800") ^ FOV of the primary focal plane covers an area of « (9mm)^ and will illuminate 43 rows x 86 columns of 
the microshutters. The rows run in the direction of the dispersion, so for unconfused spectral acquistion only one shutter 
per row can be open. 

The large number of slitlets presents a formidable target acquisition challenge, which is made more pressing given the 
short time available for target selection during a sounding rocket flight, with a total observing time ^ of 400 s. Previous 
knowledge of the bright regions in the target field could be used to pre-program the array. However, this strategy could 
easily be foiled by an Acquistion and Control System (ACS) misalignment with the telescope bore center, which experience 
has shown can be as large as 5' in pitch and yaw and > 1° in roll. A robust alternative is an automatic targeting system, 
which for the sounding rocket flight will be a simple brightest-object on a row targeting algorithm. On target, the ACS 
provides excellent hold accuracy ^1" and a command uplink system is provided to fine tune pointing misalignments in 
realtime after target acquistion at w T+85 s. 

Our baseline target acquistion plan is to launch with the microshutters open. At arrival on target we will assess, 
fine tune and lock the pointing using the FORTIS zero-order image. After a set integration period (« 30 s) the onboard 
microprocessor will generate a 43 x 86 brightest row member bitmask from the accumulated zero-order image and pass it 
to the 584 shift registers. The waveform will be synchronized with the magnet sweep, using a field programable gate array 
(FPGA) to close all but the selected shutters. Unconfused spectral observations will then commence for the duration of the 
flight, ending at « Th-530 s. 

8. CURRENT PROGRESS ON THE FORTIS BUILD 

We have purchased and are testing the 0.5 m primary mirror. The tri-axial elliptical secondary is in fabrication. A spherical 
test secondary mirror has been delivered to the grating vendor and a proper elliptical secondary has been delivered to JHU 
as a comparison against which to test the tri-axial optic. The detector is in fabrication. The optical bench is in design 
phase along with the zero-order microshutter interface. Microshutter arrays have been fabricated but await selection. We 
anticipate integration of the spectro/telescope will begin in the summer of 2009 and expect to begin integration of FORTIS 
with the NASA sounding rocket launch provider in summer of 2010. Our first target is currently scheduled to be M33. 
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